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We have studied the structure and thermodynamic properties of isotropic three-dimensional core-softened 
fluid by using the second-order Ornstein-Zernike integral equations completed by the hypernetted chain and 
Percus-Yevick closures. The radial distribution functions are compared with those from singlet integral equa- 
tions and with computer simulation data. The limits of the region of density anomaly resulting from different 
approximate theories are established. The obtained results show that the second-order hypernetted chain ap- 
proximation can be used to describe both the structure and the density anomaly of this model fluid. Moreover, 
we present the results of calculations of the bridge functions. 
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1. Introduction 

Recently, significant research activity has been focused on the so-called core-softened models 
in which the repulsive part of the intermolecular interaction potential exhibits a softening region 
at "intermediate" interparticle separations, in addition to usual hard or soft (as in the case of 
Lennard- Jones (12,6) potential) repulsive branch at short distances. The softening region can be 
described by a linear or nonlinear ramp, a shoulder, a single or multiple attractive well, or a 
combination of all these features jlUllj. Model fluids with spherically symmetric core-softened 
potentials exhibit anomalous thermodynamic and dynamic behaviours [12Hl4j that usually occur 
in real fluids with directional interparticle interactions, e.g. water, silica, some liquid metals and 
phosphorus. Examples of the anomalous behavior include the existence of a density maximum as 
a function of temperature, an increase of the diffusion coefficient upon compression, and, for some 
systems, the existence of multiple fluid-fluid transitions. 

The reason why the soft-core model potentials can reproduce the density anomaly is that at 
low densities and temperatures, the neighboring molecules are separated by distances coinciding 
with the range of attractive potential well to minimize the energy. At higher temperatures, the 
particles can penetrate into the energetically less favorable softened core to gain more entropy, 
thus giving rise to an anomalous contraction upon heating, see e.g., [la. 1611 . The systems with 
core-softened potential have been studied using computer simulations [442 1| as well as integral 
equations methods @-|<J H EH] • 



Amongdifferent models of the core-softened potentials, the potential used by Barros de Oliveira 
et al. [22h24[ | that yields both the liquid-gas phase transition and thermodynamic, structural and 
dynamic anomalies seems to be interesting. The potential is the sum of a Lennard- Jones potential 
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of the well depth Sf and a Gaussian barrier centered on radius r = ?*o with the height of ast, 
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Depending on the choice of parameters, equation (1) represents the whole family of two length 
scales intermolecular interactions, from a deep double-well potential 0, [2^, [2(| to a repulsive 
shoulder [J, 12ZL \2M- Moreover, due to its continuity, the potential (1) is convenient both for the 
Monte Carlo and molecular dynamics simulations. 

For specific choices of the parameters, the potential (1) can exhibit a double well behavior, 
similar to the potential studied by Cho et al. [3j, [25J. However, the attractive double well may bring 
both the liquid-gas phase transition and the anomalies to higher temperatures into an unstable 
region of the phase diagram in the pressure-temperature plane [26| . In order to circumvent this 
behavior, several studies were carried out for a specific set of parameters that give 

rise to a potential with a repulsive ramp followed by a small attractive well. In particular, the 
following values were often used: a = 5, c = 1 and ro = 0.7crf . It was shown [22 - 24 1 that for the 
above set of parameters, the pressure-temperature curves have a minimum if calculated at certain 
densities (i.e. (dP/dT) p — 0). Consequently, the derivative (dp/dT)p = for some thermodynamic 
states. The region of density anomaly corresponds to state points for which (dp/dT)p > and is 
bounded by the locus of points for which the thermal expansion coefficient is equal to zero. 

In our recent work |3l| we applied grand canonical Monte Carlo simulation and integral equa- 
tions with hypernetted chain (HNC), as well as Rogers- Young closures to study the above defined 
model. We found still another anomaly in the system, namely we demonstrated that for some tem- 
peratures the derivative of the density with respect to the chemical potential exhibits a minimum 
at a certain density, followed by a maximum. Also, a peculiarity in the dependence of the specific 
heat upon density was found. We compared the pair distribution functions resulting from integral 
equations and from computer simulations and established that common HNC approximation was 
not successful in capturing the region of anomalies in contrast to Rogers- Young approximation that 
imposes thermodynamic self-consistency by construction. However, the HNC was very accurate at 
high fluid densities. It is worth mentioning that a popular mean spherical approximation has not 
been used to study the model in question because this approach intrinsically requires an adequate 
splitting of the potential into a short- and long-ranged parts. This approximation does not seem 
to be beneficial for some soft-core models with two length scales and of this study in particular. 

The results of integral equations reported in the literature for the potential (1) were based on 
the singlet OZ equation (33| . A more sophisticated approach, but one that is also more demanding 
of computational resources, is based on the inhomogeneous OZ equation [34| . 

The inhomogeneous OZ equation has usually been used not only to calculate the structure 
and thermodynamic properties of fluids in contact with surfaces [35H39(, but also to determine 
the structure and the surface tension at the liquid- vapor interface [40]. However, one can also 
assume that for a single-component fluid, the source of an external potential field is just a single 
distinguished particle identical to all remaining molecules of the system. This is the so-called Percus' 
trick [41, 42 1 or the "test particle" method. The method was extended later by Attard 4^, 44 1 in 
the framework of the second-order OZ approach for bulk fluids, see e.g. [45l - l47| for the discussion of 
related important issues. Moreover, this approach was successfully applied to several model fluids 
with spherically symmetric intermolecular potentials (48- 54]. Better description of the desired 
properties or in some cases even qualitatively new findings in comparison to the singlet level 
theory were obtained. However, to our best knowledge, this approach has never been tested for 
the core-softened models of the interparticle potentials. Therefore, the motivation of this work is 
to use the second-order integral equations to the fluid of particles interacting via the core-softened 
potential (1). The second-order equations are applied to the study of the microscopic structure 
and thermodynamic properties of the fluid, in particular, in view of its anomalous behavior in a 
certain region of thermodynamic states. The results are compared with computer simulation data 
and with the predictions of the singlet integral equations [3l|. We explore here the second-order 
Percus- Yevick and hypernetted closures to the nonuniform OZ equation. 
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2. Theory 

The OZ equation for an inhomogeneous fluid, wherein the density p(r) is not constant reads 

/*2(ri,r 2 ) = C2(ri,r 2 ) + J dr 3 /i 2 (ri, r 3 )c 2 (r 3 , r 2 )p(r 3 ), (2) 

where /i 2 (ri,r 2 ) and c 2 (ri,r 2 ) are the total and direct correlation functions of an nonuniform 
fluid. We have used here the subscript "2" in order to distinguish these functions from the common 
uniform fluid correlation functions. 

Generally, equation @ applies when the inhomogeneity is due to an external potential field. 
However, one distinguished molecule can be also considered as a source of the external potential. 
In such cases equation @ can be solved using any of the common approximations, as HNC 

/&2(ri,r 2 ) = exp[/i 2 (ri,r 2 ) - c 2 (ri,r 2 ) - /3«(|ri - r 2 Q] - 1, (3) 

or Percus-Yevick (PY) approximation 

c 2 (ri,r 2 ) = {l-exp[/3u(|r 1 -r 2 |)]}[/i 2 (r 1 ,r 2 ) + l]. (4) 

In the above m(|i"i — r 2 |) is the pair potential and /? = 1/kT is the inverse temperature. However, we 
stress that the above closures do not relate uniform, but rather nonuniform correlation functions 
and therefore we refer to the results of equation © with the closure given by equation or by 
equation (HI), as to the HNC2 and PY2 results, respectively, in contrast to the usual HNC and PY 
theories [33j. 

To have a set of equations completed we also need a relation between p(r) and the pair correla- 
tion functions. The exact equation developed by Lovett, Mou, Buff and Wertheim (LMBW) f55l.l56| 
reads 

Vb%(n)]= /dr 2 c 2 (r 1 ,r 2 )Vp(r 2 ), (5) 

where y(i"i) is the one-particle background (cavity) function, p(ri) = cxp[— f3v(ri)]y(ri) and v (ri) 
is the potential due to the particle which is regarded as the source of the inhomogeneity. In our 
case, is just the pair potential, v(ri) = — 0|), where the distinguished particle is set at 

the origin, r 2 = 0. 

One should mention here an important difference between the theories used to study systems 
in contact with a surface [35M39l | (or the systems involving a gas- liquid interface and theories 
based on the method of Attard. In the former case the two-particle correlation functions reduce 
to the functions for bulk fluids, providing the both particles to be located far away from an 
inhomogeneous region. In the Attard's approach, the function /i 2 corresponds rather to the third- 
order (conditional) correlation function of a bulk fluid, while the "usual" pair correlation function 
is related to the local density. Indeed, because the bulk pair distribution function, g(r), gives the 
probability density of finding a pair of fluid molecules at a separation, r, the g(r) is related to the 
local density, p(r), via 

9(r)=p(\r\)/p, (6) 

where |r| is the distance from the distinguished particles that is the source of the system inhomo- 
geneity and p = limiri^oo p(|r|) is the bulk density. Of course, the cavity function y(r) that enters 
equation ([3]) satisfies the relation g(r) = exp[— (3u(r)]y(r). 

To solve equations @ and with the closure given by equation or equation @ we use 
the numerical algorithm Attard's [43| that relies on the expansion of the two-particle functions (ft, 2 
and c 2 ) in series of Legendre polynomials, for details see [43j. 

Note that the closure equations §3§ and ^ constitute an approximation, similar to the common 
bulk theory, whereas the OZ equation ([2]) and the equation for the profile, ([5]) are exact. 

Equation (JSJ) satisfies the relation g(r) — exp[— (3u{r)]y{r) . The bulk direct correlation function, 
c(r), can be recovered from g{r) using the bulk Ornstein-Zernike (OZ) equation 

h(\r 12 \) - cflrial) = p [ dr 3 MM)c(|r 23 |), (7) 
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where h(r) = g(r) — 1. 

The knowledge of the total, h(r), and the cavity, y(r), functions enable us to determines the 
bridge function, B(r). The bridge function is defined as, see e.g. [33, 57| 



B(r)=lny(r)- 7 (r), (8) 

where j(r) — h(r) — c(r). The bridge function plays a key role in theories based on the singlet OZ 
equation (7). It is given by the smallest set of diagrams in the graphical expansion, and all other 
functions can be calculated from it via the OZ equation [33j | . 

The main difference between the first- and the second-order integral equation theories for bulk 
fluids is that in the latter approach, the bridge function results from the calculations performed, 
whereas in the case of first-order theories the bridge function is imposed as a closure. Previous 
calculations carried out for hard spheres, as well as for Lennard- Jones fluids, gave insight into the 
course of the bridge functions. However, to our best knowledge no results of calculations of the 
bridge function for core-softened potential models have been presented so far. 

3. Results and discussion 

The second-order integral equations were solved using the Attard's algorithm [43j. Usually, 80 
Legendre polynomials were used and the grid in r was 0.04cr. To test the accuracy we increased 
the number of Legendre polynomials to 120 and decreased the grid size to 0.025ct. Moreover, we 
also solved the singlet integral equations and carried out grand canonical ensemble Monte Carlo 
simulations. The details on the two latter methods are presented in our previous work [3l| . We 
also note that the pressures given below were calculated from the virial equation of state (see 
equation (6) of [H|). 

Similarly to our previous work[3f], all the calculations have been carried out assuming that 
the parameters of the potential (1) are a = 5, c = 1 and ro = 0.7<7f . First, we concentrate 
on a comparison of the pair distribution functions obtained from the second-order OZ equation, 



the singlet OZ equation and the Grand Canonical ensemble (GCMC) computer simulations [3l|. 
Comparisons are for thermodynamic conditions that cover the density anomaly region found in 
molecular dynamics studies [23]. In figures 1 (a)-(c) we display the functions g(r) evaluated at the 
density p* = paf — 0.1 and at three temperatures, T* — kT/ef — 0.5, 0.3 and 0.2. At the highest 
considered temperature the results of the singlet theories (PY, HNC and Rogers- Young, (RY)), the 
second-order theories, PY2 and HNC2 are compared with GCMC data. Before discussing the results 
shown in figure 1, we recall that the singlet RY closure contains one adjustable parameter and this 
parameter is evaluated requiring the equality of pressures from the virial and compressibility routes, 
for details see |3lj . 

In figure 1 (a) we see that the singlet PY approximation yields the radial distribution function 
that significantly deviates from all remaining results. The second-order PY2 approximation is def- 
initely more accurate than the singlet PY theory, but its results still differ much from computer 
simulations. In contrast to the singlet PY theory, the singlet HNC closure yields the results that 
qualitatively agree with computer simulations, but among all the singlet theories the thermody- 
namically self-consistent RY closure gives the best predictions. For the thermodynamic state from 
figure 1 (a) the second-order hypernetted chain approximation reproduces the GCMC data with 
the highest precision. We stress that in contrast to the singlet RY theory the HNC2 approach 
contains no adjustable parameter. 

The accuracy of the second-order HNC2 theory is also confirmed by the results presented in 
figure I (b). HNC2 theory well reproduces the location and the height of the first maximum of 
g(r) and the shoulder at distances r/at < 1.4. However, the first minimum of g(r) at r/df « 3.32 
is a bit better predicted by the RY approximation. Similar conclusions can be drawn from the 
inspection of the curves presented in figure 1 (c). Since the PY2 approximation actually fails to 
reproduce the simulation data with a reasonable accuracy, the PY and PY2 results are omitted in 
figure 1 (c) and in the following figures 2 and 3. The singlet PY approximation is quite accurate 
when applied to a hard-sphere system [33]. On the contrary, in the case of soft potentials, e.g., 
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Figure 1. (Color on-line) Radial distribution functions evaluated from the singlet (PY, HNC 
and RY), second-order (PY2, HNC2) theories and from GCMC simulation at p* =0.1 and at 
T* — 0.5 (part a), 0.3 (part b) and 0.2 (part c). The nomenclature of the lines is given in part a. 
Note that consecutive parts b and c contain the results of selected theories. 



Gaussian- like potential, the singlet HNC approximation works quite satisfactory [58| • It is thus not 
surprising that for core-softened potential, the hypernetted chain approximation (at singlet and 
second-order levels) performs better than the the Percus-Yevick closures. However, the observed 
big differences between the predictions of these approximations are surprising, because in the case 
of the systems studied so far (IjWHsj the differences between the PY2 and HNC2 approximations 
were rather small [59|. 

Further tests of the accuracy of the HNC2 approximation are shown in figures 2 and 3. Figure 2 
shows the results at the density p* — 0.14, consecutive parts a - dare for the temperatures T* = 0.5, 
0.3, 0.2 and 0.15, respectively. Again, the HNC2 approximation is superior over the best singlet RY 
theory. In particular, it better describes the formation of the first peak of g(r) within the repulsive 
ramp of the potential (1) (the plot of the potential (1) gives figure 1 of [23]), though both RY and 
HNC2 theories underestimate the height of this maximum. By contrast, the singlet HNC approach 
describes much higher first peak of g(r) compared to GCMC data. 

Finally, we calculated the radial distribution functions at higher densities, p* = 0.16 and 0.20 
(figure 3). For these densities, the second-order HNC2 equation performs well. In particular, its 
predictions of the height of consecutive maxima and minima of g(r) are more accurate than the 
singlet RY approximation. However, for some state points, the accuracy of the singlet HNC equation 
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Figure 2. (Color on-line) Radial distribution functions evaluated from the singlet (HNC and 
RY), second -order (HNC2) theories and from GCMC simulation at p* = 0.14 and at T* = 0.5 
(part a), 0.3 (part b), 0.2 (part c) and 0.15 (part d). The nomenclature of the lines is given in 
figure 1 (a). 

is quite good and even better than the HNC2 theory, cf. figures 3 (b) and 3 (c). 

Before discussing an anomalous behavior of the system in the pressure-temperature plane, that 
results from the second-order integral equations, we briefly recall previous findings. The line of 
temperatures of maximum density (TMD) was determined from NPT molecular dynamics (23| and 
singlet PY [23| and RY [23l[3l| integral equations. It was also demonstrated (23l.l3l1| that the singlet 
HNC approach fails to predict this anomaly. 

Molecular dynamic simulations revealed that for densities 0.12 < p* < 0.14, the pressure- 
temperature curves at constant density have a minimum (see also our comment in the introductory 
section) which implies density anomaly. The maximum temperature at which the density anomaly 
still exists is Tj* ax rs 0.25 [23|. The predictions of the singlet RY approximation depend on the 
method according to which its adjustable parameter was determined (3l| . Barros de Oliveira et 
al. [2^1 obtained the adjustable parameter of the RY equation by checking the consistency between 
the compressibilities, calculated from the virial and compressibility equations of state. However, 
in our work [3l| we imposed the so-called global coexistence criterion, according to which the 
pressures from two above mentioned thermodynamic routes were compared. Barros de Oliveira et 
al. [2^ found that the density anomaly exists for the densities from the range of [0.12,0.14] and 
that T^ ax ~ 0.23. Our calculations, however, led to a wider interval of the densities, [0.12, 0.19]. We 
should also note that the singlet PY approximation predicts (23| much wider region of densities, at 
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Figure 3. (Color on-line) Radial distribution functions evaluated from the singlet (HNC and 
RY), second -order (HNC2) theories and from GCMC simulation at p* = 0.16 and at T* = 0.5 
(part a), p* — 0.16 and at T* — 0.3 (part b), p* = 0.2 and at T* = 0.5 (part c), and p* — 0.2 
and at T* = 0.3 (part d). The nomenclature of the lines is given in figure 1 (a). 



which the density anomaly occurs, [0.13, 0.3]. Instantaneously, this region is shifted towards much 
higher temperatures, with T* lax = 0.86. 

Figure 4 shows the isochores evaluated from different theories (the reduced pressure is defined as 
P* = Pafei). At low densities {p* < 0.07; the relevant curves are not shown for the sake of brevity) 
none of the approximations predicts the existence of the density anomaly. For p* = 0.08 (part a) 
only the PY2 theory yields a weak minimum of the pressure, all the remaining approximations do 
not show any anomaly. However, it is interesting to note that the results of the singlet PY and 
HNC2 approximations almost coincide, except for low temperatures. For p* =0.1 (part b) only the 
PY2 isochore exhibits anomaly, but now the minimum of the pressure is preceded by a maximum. 
This behavior is completely unexpected, but, it is rather an artifact of the theory, because no 
computer simulations [23| provide its confirmation. For p* = 0.14 (part c) all the approximations, 
but the HNC1 predict the density anomaly, whereas at higher bulk density, p* = 0.2, (part d) only 
PY1 and PY2 predict anomalies. 

The range of the densities for which the HNC2 leads to anomalous behavior is rather narrow, 
0.13 ^ p* ^ 0.155. As we have shown in our previous work [3l|, the RY approximation with the 
global coexistence criterion predicts anomalous behavior up to p* « 0.19. Surprisingly, the PY2 
approximation yields two density anomaly regions: the first one is for the densities [0.07,0.11] and 
the second one is for p* ^ 0.13. We must recall, however, that the structure predicted from PY2 
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Figure 4. (Color on-line.) Isochores resulting from different approximate theories for different 
densities, given in consecutive parts. The nomenclature of the lines is displayed in part a. 




Figure 5. (Color on-line) The TMD lines from different approximate theories and from MD 
simulations [23| (open circles). Abbreviations are as follows: HNC2 denotes the second-order 
hypernetted chain theory results; PY is the singlet Percus-Yevick theory [23^ . RY denote the 
Roger- Young results obtained by Barros de Oliveira et al. and RY-gcc are the Roger- Young 
results calculated using the global coexistence criterion. The inset magnifies a part of the main 
figure. Dotted lines here indicate the thermodynamic conditions at which the bridge functions 
displayed in figure 6 were calculated. 
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approximation differs very much from the computer simulation results. Summary of our calculations 
is presented in figure 5. We show here TMD lines resulting from different approximate theories. 
Since the PY2 predictions quantitatively differ from the computer simulation results, the relevant 
curve has been omitted. The TMD line from the RY approach evaluated using global coexistence 
criterion differs from that obtained by Barros de Oliveira et al. [23j None of the approximate 
theories is capable of reproducing the simulation results at quantitative level. 

Finally, we show the results of the bridge function calculations. The bridge function for the 
model in question has never been investigated so far. The function 7(7"), entering equation (8) was 
calculated from 



7(r) 



1 



2tt 2 ? 



dfcfc sin(fcr) 



ph(k) 2 
l + ph{k) 



(9) 



where h(k) is the Fourier transform of the function h(r). We should stress that the calculations of 
the Fourier transform h(k) and the inverse transform in equation (9) should be carried out with 
a special care, as described by Kolafa et al. [60]. Our calculations were performed at a constant 
density equal to p* = 0.152 for a set of temperatures across the TMD line, as well as at a constant 
temperature, T* — 0.175 and for several densities, again across the TMD line (cf. figure 5). The 
functions B(r) significantly differ from the functions for hard-spheres (48l [60|. In particular, the 
decay of the bridge functions for the core-softened potential is much slower than for hard-spheres. 
Even for r/crt ~ 10, the bridge functions exhibit well pronounced oscillations. These oscillations 
almost vanish at a distance as large as r/<7f w 16, whereas for hard-spheres at very high densities, 
the oscillations already vanish at r/Vf « 5, cf. figure 3 of [6(|- Crossing the TMD line seems to 
have no effect on the shape of the bridge functions. The evolution of the functions B(r) along the 
temperature and the density branches (figure 6) is smooth and we do not observe any peculiarities 
connected with specific thermodynamic behavior of the system. However, the range of oscillations 
of a bridge function is worth stressing once again. 
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Figure 6. (Color on-line) Changes of the bridge functions across the TMD line. Part a is at a 
fixed density, p* = 0.152, and for different temperatures given in the figure, whereas part b is at 
a fixed temperature, T* = 0.175, and for three different densities given in the figure. The paths 
along which the calculations were performed are shown in the inset to figure 5 as dotted lines. 

Let us briefly summarize our findings. We have demonstrated that the radial distribution func- 
tions predicted by the second-order HNC2 approach are more accurate than the predictions of any 
singlet theory. In contrast to the singlet RY approach, the HNC2 approximation does not involve 
any adjustable parameter. However, the PY2 approximation is inaccurate and the distribution func- 
tions resulting from it significantly differ from those of all remaining theories. The discrepancies 
between PY2 and HNC2 approximations are puzzling and we cannot offer any physical interpre- 
tation of these trends. The second-order HNC2 approach yields the density anomaly in the system 
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for densities from the interval [0.13, 0.15]. The maximum temperature at which this phenomenon is 
observed is T* „ « 0.2. Moreover, for the first time we have evaluated the bridge functions for the 
core softened fluid. Unlike the hard-sphere or Lennard- Jones fluids, the oscillations of the bridge 
functions extend over much larger interparticle separations. A successful parametrization of the 
bridge function similar to hard-sphere systems, see e.g. [60j-|62|, would be helpful in developing new 
closures for singlet level theory for a class of core-softened potentials. 



References 

1. Hemmer P.C., Stell G., Phys. Rev. Lett., 1970, 24, 1284; doi|10.1103/PhysRevLett.24.1284[ 



2. Stillinger F.H., Head-Gordon T., Phys. Rev. E, 1993, 47, 2484; doi: 10.1103/PhysRevE.47.2484 

3. Cho C.H., Singh S., Robinson G .W., Phys. Rev. Lett., 1996, 76, 1651; 
doi jlO. 1 103 /PhysRevLett .76. 1651 1 

4. Jagla E.A., Phys. Rev. E, 1998, 58, 1478; doi |10.1103/PhysRevE.58.1478| 

5. Sa dr-Lahijany M.R., Scala A., Buldyrev S.V., Stanley H.E., Phys. Rev. Lett., 1998, 81, 4895; 
doi:| 10. 1 103 /PhysRevLett .81 .4895 1 



6. Egorov S.A., J. Chem. Phys., 2008, 129, 024514; doi: 10.1063/1.2953332 



7. Malescio G., Pellicane G., Phys. Rev. E, 2001, 63, 020501(R); doi: 10.1103/PhysRevE.63.020501 

8. Kumar P. et al., Phys. Rev. E, 2005, 72, 021501; doi: 10.1103/PhysRevE.72.021501 



9. Chakraborty S.N., Chakravarty C, Phys. Rev. E, 2007, 76, 011201; doi: 10. 1103/PhysRevE.76. 011201 

10. Zhou S., Lajovic A., Jamnik A., J. Chem. Phys., 2008, 129, 124503; doi |10.1063/1.2982162 

11. Chaimovich A., Shell M.S., Phys. Chem. Chem. Phys., 2009, 11, 1901; doi: 10.1039/b818512c 



12. Wilding N.B., Magee J.E., Phys. Rev. E, 2002, 66, 031509; doi: 10.1103/PhysRevE.66. 031509 

13. Gibson H.M., Wilding N.B., Phys. Rev. E, 2006, 73 061507; doi|10.1103/PhysRevE.74.019903 

14. Yan Z.Y. et al, Phys. Rev. E, 2006, 73, 051204; doi |10.1103/PhysRevE.73.051204| 

15. Krekelberg W.P., Mittal J., Ganesan V., Truskett T.M., Phys. Rev. E, 2008, 77 041201; 



62[~ 
512cl 
1509] 
119903 



doi: 10.1103/PhysRevE.77.041201 



16. Ba rros de Oliveira A., Netz PA. , Barbosa M.C., Eur. Phys. J. B, 2008, 64, 481; 
doi (iail40/epjb/e 2008-0010T6] 

17. Gribova N.V., Fomin Yu.D., Frenkel D., Ryzhov V.N., Phys. Rev. E, 2009, 79, 051202; 



doi: 10. 1103/PhysRevE.79. 051202 



18. Lomba E., Almarza N.G., Martin C, McBride C, J. Chem. Phys., 2007, 126, 244510; 
doi |10. 1063/1.2748043| 

19. Franzese G. et al., Phys. Rev. E, 2002, 66, 051206; doi |10.1103/PhysRevE.66.051206 

20. Malescio G. et al., Phys. Rev. E, 2005, 71, 061504; doi |l0.1103/PhysRevE.71.061504 

21. Barraz N.M., Salcedo E., Barbosa M.C., J. Chem. Phys., 2009, 131, 094504; doi |10.1063/1.3213615 

22. Ba rros de Oliveira A Netz PA., Colla T., Barbosa M.C., J. Chem. Phys., 2006, 125, 124503 
doi |i"ai063/l .2357119] 

23. Ba rros de Oliveira A., Netz PA., Colla T., Barbosa M.C., J. Chem. Phys., 2006, 124, 084505 
doi fiai 063/1. 2168458] 

24. Ba rros de Oliveira A., Barbosa M .C., Netz PA., Physica A, 2007, 386, 744; 
doi|10.1016/j.physa.2007.07.0T5~l 



25. Cho C.H., Singh S., Robinson G.W., J. Chem. Phys., 1997, 107, 7979; doi 10.1063/1.475060 

26. Ne tz P.A., Raymundi J.F., Cam era A.S., Barbosa M.C., Physica A, 2004, 342, 48; 
doi |10.1016/j.physa.2004.04.058| 

27. Jagla E.A., J. Chem. Phys., 1999, 110, 451; 1999, 111, 8980; doi |10.1063/1.478105| 

28. Jagla E.A., Phys. Rev. E, 2001, 63, 061501; 2001, 63, 061509; doi |10.1103/PhysRevE.63.061501| 

29. Dominguez H., Pizio O., Pusztai L., Sokolowski S., Adsorpt. Sci. Technol., 2007, 25, 479. 

30. Piz io P., Dominguez H., Puszta i L., Sokolowski S., Physica A, 2009, 388, 2278; 
doi |iai016/j.physa.2009.03.014l 

31. Piz io P., Dominguez H ., Duda Yu., Sokolowski S., J. Chem. Phys., 2009, 130 174504; 
doi flOl 063/1. 3125930] 

32. Ba rros de Pliveira A., Salcedo E., Chakravarty C, Barbosa M.C., J. Chem. Phys., 2010, 132, 234509; 
doi (l0l 063/1 .3429254] 

33. Hansen J. P., McDonald I.R., Theory of Simple Liquids. Academic Press, 3rd edition, New York, 2006. 

34. Henderson D., Sokolowski S., Borowko M., in: Computational Mathods in Surface and Colloid Science, 
ed. M. Borowko. Marcel Dekker, New York, 2000. 



13601-10 



Core-softened model fluid 



35. Sokolowski S., J. Chem. Phys., 1980, 73, 3507; doi: 10.1063/1.440507 

36. Sokolowski S., Mol. Phys, 1983, 49, 1481; doi: 10.1080/002689783001 02091] 

37. Plischke M., Henderson D., J. Phys. Chem., 1984, 88, 6544; doi |10.1021/jl50670a015"1 

38. Henderson D., Sokolowski S., Trokhymchuk A., J. Chem. Phys., 1995, 103, 4693; doi |10.1063/1.470656 

39. He nderson P., Pizio P., Trokhymch uk A., Sokolowski S., Physica A, 1997, 244, 147; 
doi|10.1016/S0378-4371(97)00237-9| 



40. Omelyan I., Hirata F., Kovalenko A., Chem. Phys. Phys. Chem., 2005, 7, 4132; doi: 10.1039/b507761c 

41. Percus J.K., Phys. Rev. Lett., 1962, 8, 462; doi |10.1103/PhysRevLett.8.462| 

42. Percus J.K. - In: The Equlibrium Theory of Classical Fluids, Eds. H.L. Frisch and J.L. Lebowitz. 
Benjamin, New York, 1964, 33. 

43. Attard P., J. Chem. Phys., 1989, 91, 3072; doi| 10.1063/l~4 56930 

44. Attard P., J. Chem. Phys., 1989, 91, 3083; doi |l0.1063/1.45693l| 

45. Henderson D., Condens. Matter Phys., 1993, 2, 31. 

46. Holovko M.F., Pizio O. Preprint ICMP-91-2U, 1991 (in Ukrainian). 

47. Holovko M.F., Pizio O., Haymet A.D.J. , Henderson D. Preprint IT P-91-55E, 1991. 

48. Henderson D., Sokolowski S., J. Chem. Phys., 1995, 103, 7541; doij lO. 1063/1.470322 

49. Henderson D., Sokolowski S., J. Chem. Phys., 1996, 104, 2971; doi |l0.1063/1.471118 

50. Henderson D., Sokolowski S., Condens. Matter Phys., 1996, 7, 51. 

51. Du h P.M., Henderson P., Mier y Teran L., Sokolowski S., Mol. Phys., 1997, 90, 563; 
doi |10. 1080 /0026897970 9482638 

52. Henderson P., Sokolowski S., Mol. Phys., 1997, 90, 85; doi | 10. 1080 /002689797172895 1 

53. Henderson P., Sokolowski S., Wasan P., Mol. Phys., 1998, 93, 294; doi |10.1080 /002689798169302 

54. Bryk P., Pizio O., Sokolowski S., J. Chem. Soc. Faraday Trans., 199 7, 93, 2367; doi|10 1039/a608565b"1 

55. Lovett R., Mou C.Y., Buff F. P., J. Chem. Phys., 19 76, 65, 570; doijlO 1063/1.433110[ ~ 

56. Wertheim M.S., J. Chem. Phys., 1976, 65, 2377; doi |10.1063/1.433352| 

57. Yukhnovski I.R., Holovko M.F., Statistical Theory of Classical Equilibrium Systems. Naukova Eumka, 
Kiev, 1980. 

58. Louis A. A., Bolhuis P.G., Hansen J. P., Phys. Rev. E, 2000, 62, 7961; doi |10.1103/PhysRevE.62.7961| 

59. Henderson P. - In: Fundamentals of Inhomogeneous Fluids, ed. P. Henderson. Marcel Pekker, New 
York, 1992, Chapter 4. 

60. Kolafa J., Labfk S., Malijevsky A., Mol. Phys., 2002, 100, 26 29; doi[Tol 080/002689702 1013635"7| 

61. Grundke E.W., Henderson P., Mol. Phys., 1972, 24, 269; doi |10. 108 0/0026897720 01014311 
Henderson P., Grundke E.W., J. Chem. Phys., 1975, 63, 601 ; doi|10.1063/1.431378"T~ 

62. Labfk S., Malijevsky A., Mol. Phys., 1984, 53, 381; doi jl0.1080/00268978400102381 | 
Malijevsky A., Labfk S., Mol. Phys., 1987, 60, 663; doi: 10. 1080 /00268 978700100441 [ 



13601-11 



O. Pizio, Z. Sokotowska, S. Sokotowski 



MiKpocKoniHHa CTpyKTypa i TepMOflMHaiviiKa MOfle/ii n/iMHy 3 
noM'flKiueHMM KopoM 3 TeopiV iHTerpa/ibHMX piBi-mHb flpyroro 
nopjiflKy 

0. ni3ic0, 3. CoKO/IOBCbKfl C. CoKO/IOBCbKI/ltP 
iHCTMTyT xin/iiT yrAM, KowoKaH, MeKCWKa 

I HCTMTyT arpocpi3HKH riojibCbKoT aKafleMi'f HayK, T\\o6mH, PecnyS/iiKa l~lo.nbma 
yHiBepcuTeT iM. MapiT CK/iaflOBCbKoT-Kiopi, T\K>6n\n, Pecny6/iiKa l~lo.nbLU,a 

Mm flocnifln/in CTpyKTypHi i TepMOfli/iHaMinm B/iacTH Boc~ri oflHopiflHoro TpwBWMipHoro njiMHy 3 noM'smnieHi/iM 
KopoM, BUKopwcTOByioHW iHTerpa.nbHi piBHHHHH OpH WTefi Ha-L|epHiKe flpyroro nopaflKy i3 rinepjianmioroBi/iM 
3a m m Ka hhsm Ta 3a m n Ka h h?i m l~lepKyca-€BiKa. 3po6/ieHO nopiBHnHHfl paflia/ibHux cpyHKLiiw po3nofli.ny 3 Biflno- 

BiflHUMM CpyHKLliaMU, OTpMMaHWMM i3 CH H T/ieTH MX iHTerpa/lbHMX pi BH51 H b, a T3 KO>K 3 flaHMMW KOM n ' KDTepHOTO 

MOflejiKDBaHHa. 3 pi3Hwx Ha6jiw>KeHnx TeopiW BCTa HOB/ieHO rpam/mi o6.nacTi aHOMa/iiT rycTWHW. OTpi/iMam pe- 
3y/ibTaTW noKa3yK)Tb, mo rinep/iaHi4ioroBe Ha6/in>t<eHHa iHTerpaiibHwx piBH^Hb flpyroro nopflflKy MO>Ke 6y~n/i 
bm KopncTa He fl/ia onucy sk CTpyKTy pi/i , TaK i aHOMa/ii'f rycTHHH L4boro MOfle/ibHoro rmwHy. Kpiwi Toro, mi/i 
npeflCTaBJiaeMO pe3y/ibTaTH o6Hi/icneHb mIctkobi/ix cpyHKL4i£n. 

KyiKJHOBi c/ioBa: reop'ia piflMH, mrerpajibm p'ibhhhhh flpyroro nopaflxy, aHOMa/iix rycriAHW, mictkobi 
cpyHKijiY 
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